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Abstract. Recent detections of GeV photons in a few GRBs by Fermi-LAT have led to strong constraints on the bulk Lorentz 
factor in GRB outflows. To avoid a large yy optical depth, minimum values of the Lorentz factor are estimated to be as high as 
800-1200 in some bursts. Here we present a detailed calculation of the yy optical depth taking into account both the geometry 
and the dynamics of the jet. In the framework of the internal shock model, we compute lightcurves in different energy bands 
and the corresponding spectrum and we show how the limits on the Lorentz factor can be significantly lowered compared to 
previous estimates. 
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INTRODUCTION 

The compactness problem. The short time scales observed in GRBs (down to a few ms) can be used to deduce an 
upper limit on the size of the emitting region producing y-rays. This information combined with the huge isotropic 
y-ray luminosities deduced from the measured redshifts imply huge photon densities. Then the simplest assumption 
of an emission produced by a plasma radiating isotropically with no macroscopic motion predicts that y-ray photons 
should not escape due to yy annihilation yy — > e + e~ . This is in contradiction with the observed GRB spectra which 
are non-thermal and extend well above the rest-mass electron energy m e c 2 ~ 5 1 1 keV. Observation and theory can 
be reconciled by assuming that the emitting material is moving at ultra-relativistic velocities [1]. This is mainly due 
to relativistic beaming. First the relativistic beaming implies that the observer will see only a small fraction of the 
emitting region: the constraint on the size of the emitting region is now less severe. Second the collimation of photons 
in the same direction reduce the number of potential interactions. Finally the typical yy interaction angle becoming 
small the photon energy threshold for pair production becomes higher. This theoretical context combined with the 
observational data gives the possibility to estimate a minimum Lorentz factor r m ; n for the emitting outflow in GRBs 
[2] (or directly a Lorentz factor estimate if the yy cutoff is clearly identified in the spectrum, see [3]). 

Severe constraints on the Lorentz factor from Fermi-LAT observations. Since the launch of Fermi in June 
2008, the LAT instrument has detected high energy photons above 10 GeV in a few GRBs. The observed y-ray 
spectrum often remains consistent with a Band function covering the GBM and LAT spectral ranges without any 
evidence of a high energy cutoff which could be identified as a signature of yy e + e~. This extension by Fermi of 
the observed spectral range upper bound from 10 MeV (e.g. BATSE) to 10 GeV implies constraints on r m j n which are 
much more severe than the ones obtained previously. In a few cases r m ; n has been estimated to be of the order of 1000 
(for example: GRB 080916C - r min = 887 [4], GRB 090510 - r min = 1200 [5]). These extreme values put severe 
constraints on the physics of the central engine which should be able to strongly limit the baryon load in the outflow. 
However these r m i n values were obtained from a simplified "single zone" model where the spatial and temporal 
dependencies are averaged out. The motivation of this work is to develop a detailed approach taking into account a 
more realistic treatment of the dynamics. 
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FIGURE 1. Opacity in a single pulse - comparison with the semi-analytical work by [6]. Left panel: our model. Right panel: 
model presented in [6]. As illustrated here, our numerical model can reproduce the normalized light-curves shown in [6] (the 
observer time T and observed fluxes eF e /Fq have normalized values). The discrepancy for T < 10~ 2 is due to numerical resolution 
limitations (this corresponds to a "true" observer time t i, s < 10~ 4 s). 



COMPUTING THE 77 OPTICAL DEPTH 

General 77 opacity formula. The 77 opacity (T rr ) is given by: 

T n (E GeV ) = / dl dQ. dE n a (E)o yy (E ,y){\ - cosy) (1) 

Jle J JE c (E GeV ,\fl) 

All the physical quantities are measured in the laboratory (or source) frame. E^eV is the energy of the photon for which 
Xyy is calculated whereas E is the energy of the interacting field photon, y represents the interaction angle between 
the GeV photon and the interacting photon and a y y is the 77 interaction cross-section between these two photons. 
E c = 2(m e c 2 ) 2 j [E(j e v (1 — cosy)] is the energy threshold of the field photon above which 77 annihilation can happen. 
Finally «n is the photon field distribution [ph • cm~ 3 • erg -1 • sr~'] at a given location and time. 

The equation (1) is made of a triple integral : the dl -integration is done over the path of the GeV photon from its 
emission location to the observer, the t/i2-integration is done over the solid angle distribution of the interacting photon 
field surrounding the GeV photon whereas the liTi-integration is done over its energy distribution. The equation (1) is 
general and can be applied to any photon emitted at a given location and time with a given propagation direction in 
the GRB outflow. 



Validation of the model. The kernel of our study is the calculation of the 77 opacity created by a spherical flash. It 
is then possible to model the case of a propagating radiating spherical front (representing for example a shock wave) 
by the succession of many spherical flashes. One of the critical step is the exact calculation of the photon density «n 
taking into account all the relativistic effects. Before dealing with more complex dynamical configurations within the 
internal shock framework, the validity of our approach was tested on a simple single-pulse case with a comparison to 
the previous semi-analytic study of [6] (see fig. 1). 

APPLICATION TO INTERNAL SHOCKS 

Now the model is applied to dynamical evolutions expected in the internal shock framework, where the whole 
prompt 7-ray emission is produced by electrons accelerated by shock waves propagating within a relativistic variable 
outflow. We model the dynamics via a multiple shell model where the successive collisions between shells mimic the 
propagation of shock waves [7]. Each collision produces an elementary spherical flash: the simulated light curves 
are the result of the sum of all flashes. For each high energy photon, the 77 opacity is computed by integrating 
equation (1) from its emission location to the observer taking into account the exact radiation field «n produced by 
all the collisions in the outflow. A previous study of the 77 opacity in internal shock was made by [8]. However the 
prescription used to compute Xyy was still approximate, using the local physical conditions of the outflow where the 
high energy photon is emitted and applying them to an average formula of Xyy (as can be found in [2, 4, 5]). 




FIGURE 2. temporal and spectral evolution of the yy opacity - Example of 2 synthetic GRBs. The upper part of the three 
figures illustrates the case of a mono-pulse GRB with no additional temporal substructures whereas the bottom part show the case 
of a complex multiple-pulse GRB. In the first case the Lorentz factor distribution rises monotonously from 200 to 700, leading to 
a collision between the slow and the fast parts of the outflow. The second case is similar with the addition of a short time-scale 
variability. The kinetic power of the outflow is assumed to be constant. Left panel: the y-ray light curve as a function of the observer 
time t (, s . Middle panel: the evolution of the yy opacity seen by 10 GeV photons emitted on the line of sight as a function of the 
time t a i, s at which they are received by the observer. Right panel: y-ray spectrum integrated over the GRB emission including (solid 
line) or not including (dashed line) the yy annihilation process. 




FIGURE 3. Minimum Lorentz factor for GRB 080916C Left panel: evolution of Xyy against the mean Lorentz factor of the 
outflow r following our detailed modeling (solid line) and using an average formula [4] (dashed line). Right panel: time-integrated 
spectra obtained for different mean Lorentz factor of the outflow (the shape of the initial Lorentz factor distribution stays the same) 
and the reference spectrum without yy annihilation (dashed line). 



Effects of temporal and spectral evolution. Due to the temporal evolution of Xyy, the opacity cutoff in a time- 
integrated spectrum will be smoother than a sharp exponential decay: the cutoff transition will look more like a 
power-law steepening. The detailed model presented here is appropriate to characterize from an observer point of 
view this time evolution effect. The smoother Tyy transition is due to the fact that the time integrated spectrum is a 
superposition of instant spectra which can have different yy cutoff photon energies. This time evolution takes place 
within a same y-ray pulse, and can be even stronger in a complex burst where the light curve is made of several pulses. 
Also notice that in this case cross-interactions between different pulses can become important and strongly influence 
the time evolution of T rr : these potential cross-interactions are fully taken into account by the present model. To 
illustrate this aspect two examples of synthetic GRBs are shown in figure 2. In the case of the "mono-pulse" GRB the 
yy opacity seen by the GRB photons remains within a decade over the whole emission (except at the very beginning 
when the interacting photon field progressively builds up, see also [6]): the yy cutoff in the time-integrated spectrum 
remains close to an exponential cutoff. In the case of the "complex" multiple-pulse GRB T rr covers a range of 4 
decades and the yy cutoff is close to a power-law steepening. 



Minimum Lorentz factor: detailed modeling vs. simple estimates. One natural application of our model is the 
estimate of the minimum bulk Lorentz factor r m! „ in GRB outflows. Our model reproduces the dependence of Xyy 

on the mean Lorentz factor Y of the outflow and the energy of the high energy photon EceV'- tyy 06 ^ 2 ^ ^E G ^y + ®\ 
However the normalization we obtain is different compared to a simple "single zone" model (see fig. 3, left panel) and 
gives values of r m ,„ which are lower by a factor of 2-3 (the exact factor depending on the details of the dynamical 
features of the GRB). To illustrate this aspect with an example our approach was applied to the case of GRB 080916C. 
Using our numerical model a synthetic GRB was generated which reproduces the main observational characteristics of 
the GRB: the total radiated isotropic y-ray energy (Ei SO = 8.8 x 10 54 ergs), the spectral properties (E p , a, j3 parameters 
of the Band function), the envelop of the light curve and a short time-scale variability of 2s in the observer frame. The 
study is focused on the most constraining time bin (time bin b), during which the highest observed photon energy was 
3 GeV (16 GeV in the source rest frame). r m(n is obtained by requiring that for this photon Xyy < 1 (see figure 3, right 
panel). Following our approach we find a minimum Lorentz factor r m ,„ = 360, instead of 887 which was obtained 
from an approximate "single zone" model [4]. 



CONCLUSIONS 

The detailed 77 opacity calculation model presented in these proceedings is appropriate and accurate to study many 
aspects and consequences of 77 annihilation in GRBs. In the present work we focus on the internal shock model and 
consider the consequences and signatures that 77 opacity could have in GRB observations. The model was validated 
by comparing our results to a previous semi-analytical study [6]. 

• It is shown how a detailed calculation can predict minimum Lorentz factors r m(n which are lower by a factor of 
2-3 compared to a simplified "single zone" model where spatial and temporal dependencies are averaged out. 

• The 77 cutoff transition is also characterized in time-integrated spectra. It is usually closer to a power-law 
steepening than to a sharp exponential cutoff. The exact shape of the transition strongly depends on the details of 
the GRB dynamics. 

Other effects (Hascoet et al. in prep.) can be studied with our model: 

• The temporal evolution of Xyy during a burst could favor a delay between the MeV and GeV light curves. 

• For complex GRBs, the 77 opacity could suppress the shortest time-scale features in high energy light curves 
(above 100 MeV). 

• If MeV and GeV photons are not produced at the same location (see e.g. [9, 10]), Xyy could be further lowered, 
reducing even more the constraint on r mi „. 
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